On the extensivity of the entropy S q for N < 3 specially correlated binary subsystems 
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Many natural and artificial systems whose range of interaction is long enough are known to 
exhibit (quasi)stationary states that defy the standard, Boltzmann-Gibbs statistical mechanical 
prescriptions. For handling such anomalous systems (or at least some classes of them), nonextensive 
statistical mechanics has been proposed based on the entropy S q = k (1 — '%2™L 1 p i 9 )/(q — l), with Si = 
— kTi^ 1 pi\npi (Boltzmann-Gibbs entropy). Special collective correlations can be mathematically 
constructed such that the strictly additive entropy is now S q for an adequate value of q 7^ 1, whereas 
Boltzmann-Gibbs entropy is nonadditive. Since important classes of systems exist for which the 
strict additivity of Boltzmann-Gibbs entropy is replaced by asymptotic additivity (i.e., extensivity), 
a variety of classes are expected to exist for which the strict additivity of S q (q 7^ 1) is similarly 
replaced by asymptotic additivity (i.e., extensivity). All probabilistically well defined systems whose 

adequate entropy is Si are called extensive (or normal). They correspond to a number W e $ of 
effectively occupied states which grows exponentially with the number N of elements (or subsystems) . 
Those whose adequate entropy is S q (q 7^ 1) are currently called nonextensive (or anomalous). They 

correspond to W e ff growing like a power of N. To illustrate this scenario, recently addressed Q], 
we provide in this paper details about systems composed by N = 2, 3 two-state subsystems. 
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I. INTRODUCTION 

It is well known that, if we have a system composed by 
N statistically independent subsystems (i.e., such that all 
joint probabilities factorize into the marginal ones corre- 
sponding to each subsystem), the Boltzmann-Gibbs (BG) 
entropy Sbg = — k J2i=i Pi m P« is strictly additive, i.e., 
Sbg{A x + A 2 + ... + A N ) = S B a(Ai) + S B g{A 2 ) + ... + 
Sbg(An). A plethora of physical systems is known for 
which this remarkable property still holds asymptotically 
(N — > 00). Such is the case, for instance, of virtually all 
many-body Hamiltonian systems involving short-range 
two-body interactions. This property is called exten- 
sivity, adopting the language of thermodynamics, where 
it plays an important role. Many natural and artifi- 
cial systems exist however that do not belong to this 
class, such as many-body Hamiltonian systems involv- 
ing two-body interactions whose range of interaction is 
long enough (Newtonian gravitation is one famous exam- 
ple). Such systems are known to exhibit stationary (or 
quasi- stationary or metastable) states that defy the usual, 
BG statistical mechanical prescriptions. For handling at 
least some of such anomalous systems, a generalization 
of BG statistical mechanics, has been proposed in 1988 
, which is usually referred to as nonextensive statistical 
mechanics (see [3] for reviews). It is based on the entropy 



fc(l- 



YZiPi)l(l- 1 ) Si = Sbg), which 



generalizes the BG one. It has been shown recently that 



special collective correlations can be mathematically con- 
structed such that the entropy which is strictly additive 
is now S q for an adequate value of q ^ 1 (directly deter- 
mined by the type of correlations), whereas Sbg is non- 
additive. It is easy to imagine that, in the same way that 
important classes of systems exist for which the strict ad- 
ditivity of Sbg is replaced by just asymptotic additivity 
(i.e., extensivity), a variety of classes must exist for which 
the strict additivity of S q (g ^ 1) is similarly replaced 
by asymptotic additivity (i.e., extensivity). Such systems 
would be the object of the so-called nonextensive statisti- 
cal mechanics. Then, as a kind of bizarre linguistic twist, 
it turns out that the appropriate entropy S q for such, so- 
called nonextensive systems, is in fact expected to be 
extensive. The generic scenario is therefore as follows: 
all probabilistically well defined systems are expected to 
have an entropy which is extensive; those whose appro- 
priate entropy is Sbg ( or its associated forms, such as 
those adapted to fermions and bosons) are called exten- 
sive, and those whose appropriate entropy is S q (q ^ 1) 
(or even some other entropic form) are called nonexten- 
sive. 

A quantity X(A) associated with a system A is said 
additive (see [l| , which we closely follow here) with regard 
to a specific composition of A and B if it satisfies 



X(A + B) =X(A)+X(B) 



(1) 
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where + inside the argument of X precisely indicates that 
composition. 

If, instead of two subsystems A and B, we have N of 



2 



them (Ai, A 2 , An), then we have that 



N 



N 



xC£A i ) = ^X(A i ). 



(2) 



i=l 



i=l 



If the subsystems (e.g., just the elements of the full sys- 
tem) happen to be all equal (a quite common case), then 
we have that 



X(N) = NX(1) , 



(3) 



with the notations X(N) = X(J2? =1 A) and X(l) = 
X{A X ). 

An intimately related concept is that of extensivity. 
It appears frequently in thermodynamics and elsewhere, 
and corresponds to a weaker demand, namely that of 



hm - — ; — - < oo . 



(4) 



Clearly, all quantities that are additive with regard to a 
given composition law, also are extensive with regard to 
that same composition (and liniAr^oo X(N)/N — X(l)), 
whereas the opposite is not necessarily true. Let us apply 
these remarks to entropy. 

Boltzmann-Gibbs (BG) statistical mechanics is based 
on the entropy 



w 



Sbg = -k^Pi lnp 4 , 



(5) 



with 



w 



J2Pi = 1 ' 



i=l 



where Pi is the probability associated with the i th micro- 
scopic state of the system, and k is Boltzmann constant. 
In the particular case of equiprobability, i.e., pi = l/W 
(Vi), Eq. (5) yields the well known Boltzmann principle: 



S BG = klnW . 



(7) 



From now on, and without loss of generality, we shall 
take k equal to unity. 

Nonextensive statistical mechanics is based on the so- 
called "nonextensive" entropy S q defined as follows: 



1 v^w 



(q e R; Si = Sbg) ■ (8) 



(Later on we come back onto the denomination "nonex- 
tensive"). 

For equiprobability (i.e., pi = 1/W, Vi), Eq. (8) yields 



S q = ln q W , 
with the q-logarithm function defined as 
z 1 ^ - 1 



(9) 



ln z = 



(z e M; z > 0; lni z = lnz) . (10) 



The inverse function, the q- exponential, is given by 

e z q = [l + (l-q)z] 1 ^ 1 -^ (e z 1= e z ) (11) 

if the argument 1 + (1 — q)z is positive, and vanishes 
otherwise. Following a common usage, we shall from now 
on cease distinguishing between additive and extensive, 
and use exclusively the word extensive in the sense of 
either strictly or only asymptotically additive. 

II. N SUBSYSTEMS 
A. General considerations 

Consider a system composed by N subsystems 
A\, A2, An having respectively Wa x , Wa 2 > • • • > Wa n 
possible microstates (we only address here the basic case 
of discrete microstates). The total number of possible 
microstates for the system A\ + A2 + ... + An is then in 
principle W = W Ai +a 2 +...+a n = W Ai Wa 2 --W An - We 
emphasized the expression "in principle" because we shall 
see that a more or less severe reduction of the full phase 
space might occur in the presence of a special type of 
strong correlations between the subsystems. 

We shall use the notation p A }+ A ^+-+ A « = 
1, 2, W Al ; ii = 1, 2, Wa 2 'i •••) for the joint proba- 
bilities, hence 



W Al Wa 2 Wa n 

A 1 +A 2 + ...+A„ 

/ y / J "• / J ri\l 2 ...%N 
»1 = 1 »2 = 1 »JV=1 



1 . 



(12) 



(6) The Ai — marginal probabilities arc defined as follows: 



W A9 Wa. 



E A 1 +A 2 + ...+A N 
J J i 1 i 2 ...i N 



(13) 



hence 



=1 ijv=l 



W Al 

»i=i 



(14) 



Analogously are defined all the other N — 1 one- 
subsystem— marginal probabilities. The AiA 2 — marginal 
probabilities are defined as follows: 



At+A 2 



hence 



Wa 3 W An 

Er a 1 +a 2 +...+a 1 < 
■■■ / j Pixi 2 ...i N 



W Al Wa 2 

»i=i 12=1 



(15) 



(16) 



Similarly are defined all the other [N(N — l)/2] — 1 two- 
subsystem— marginal probabilities, as well as all the other 
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TABLE II: Joint probabilities for two binary subsystems 



TABLE I: Joint probabilities for two subsystems 



many- subsystem— marginal probabilities. The most gen- 
eral N = 2 case is indicated in Table I. 

The central point that the present paper addresses is 
whether it is or not possible to satisfy all this specific 
structure of joint and marginal probabilities, and simul- 
taneously impose the condition 

S q (A 1 + A 2 + ... + A N ) = S q (A 1 ) + S q (A 2 ) + ... + S q (A N ), 

(17) 

where S q (A 1 + A 2 + ... + An) is calculated with the 
joint probabilities, and S q (A r ) is calculated with the 
A r — marginal probabilities (r — 1,2, ...,N). To simul- 
taneously satisfy, in fact, not only Eq. (17) but also 



S q (Ai + A 2 + ... + A N ) 
[S q (A 1 +A 2 )+S q (A 1 + A 3 ) 
+ S q (A N ^ + A N )]/(N-l), 



(18) 



and all similar ones calculated with all possible combi- 
nations of many-subsystem— marginal probabilities. It is 
kind of trivial that at least one solution exists, namely 
that of mutually independent subsystems with the en- 
tropy Sbg- Indeed, we can verify that the hypothesis 



Ai+A 2 + .. 
Pi\ii...i N 



N 

n 

r=l 



Pil 



(V(i ls i 2 ,. ..,»*)) (19) 



implies Eqs. (17), (18) and all the similar ones. But it 
is by no means trivial that a different choice is possible 
involving collective correlations and a value of q different 
from unity. Paper Q answered affirmatively precisely to 
this question. We shall present here details of such solu- 
tions, namely for TV = 2 and N — 3 binary subsystems. 



B. N — 2 specially correlated binary systems 

The most general joint probabilities for two binary sub- 
systems (noted A and B, with Wa = Wb = 2) are indi- 
cated in Table II. 

It trivially verifies Eq. (17) with q = 1 in Table III 
(top) and it was shown in [l| that the Table III (middle) 



possible to interpolate between Table III (top) and Table 
III (middle) through Table III (bottom) (0 < q < 1), 
where the function f q (x,y) is defined as follows: 

( P fr + (i-pf) q + (p?r + (i- P ?) q 
-[f q (ptpi)\ q 
f q (p?,pf)] q -\pf-f q (j>t,p? w 



[pi 



We verify that 



pi~P? + f q ( P i,P?)] q = i 



(20) 



fM,p?) = f q (pf,pf), 

fq(P,l) = P, 

h( P i,p?) - pip?, 

fo( P i,P?) 



Pi 



pf-i, 



(21) 



and also Eq. (17). Typical examples of the function 
f q (x,y) are shown in Fig. 1. 

Let us consider the simple case where A — B (hence 
Pi = Pi = p) ■ Tables III become respectively Tables IV, 
where f q (p) = f q (p,p) satisfies the relation 

2p«+2(l-p)«-(/,)«-2(p-/ g )«-(l-2p+/ g )« = 1. (22) 

In Fig. 2 we present f q (p) for typical values of q and 
the q— dependence of f q (l/2). It can be straighforwardly 
verified that S q (2) = 2S q (l) (see Fig. 3). 

It is instructive to see the joint probabilities of this 
simple case normalized by those of the independent case. 
The results are shown in Figs. 4. 



(top) 
(middle) 
(bottom) 



J.1) 
Pn 
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P~ fqip) 



p(l -p) 
(1-pf 



p(l-p) p(l 

l-2 P + f q (p) 
(l-pf ■ 



P) 



(23) 
(24) 
(25) 



satisfies Eq. (17) with q = with pf 



P? 



> 1. It 



We verify that q decreasing from 1 to zero inhibits the 
occupation of the states (11) and (22), and consistently 
enhances the occupation of the states (12) and (21). 



C. TV = 3 specially correlated binary systems 

The most general joint probabilities for three binary 
subsystems (noted A, B and C, with Wa = Wb = Wc = 
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TABLE III: Joint probabilities for independent two binary 
subsystems (top), for specially correlated two binary subsys- 
tems (middle), and for an interpolation between them (bot- 
tom). 
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FIG. 1: Typical surfaces for N — 2 case: q = (top), q — 0.5 
(middle; the graphically unperfect match at pf = and at 
pf = pf = 1 is only due to numerical imprecision) , and q — 1 
(bottom). 




TABLE IV: Joint probabilities for independent two binary 
subsystems (top) , for specially correlated two binary subsys- 
tems (middle), and for an interpolation between them (bot- 
tom). For each cases, pi = pf = p. 



FIG. 2: N = 2 : A = B case with Wa = W B = 2. left: 
Functions f q (p) = Pn'(p) for typical values of g. Right: in- 
dependence of / g (1/2) (/i(l/2) = 1/4). 



FIG. 3: N — 2 . p— dependences of 2S q (l) (solid curve), 
and of S q (2) (dots) for typical values of q. For q = there 
is a cutoff probability p* = 1/2 indicated by an arrow. Since 
numerical calculation of S q (2) becomes unstable when p is too 
small, we only show values of p whose corresponding entropies 
are reliably calculated. For q — 0, the rest of the curve (i.e., 
for < p < 1/2) can be obtained by using the same solution 
but permutating p and 1 — p. 

2) are indicated in Table V where the numbers without 
parentheses correspond to system C in state 1, and the 
numbers within parentheses correspond to system C in 
state 2. 
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TABLE V: Joint probabilities for three binary subsystems 
(the quantities between parentheses correspond to C being in 
state 2, whereas the others correspond to C being in state 1). 

The corresponding AB— marginal probabilities arc in- 
dicated in Table VI which precisely reproduces the situ- 
ation we had for the two-subsystem (^4 + B) problem. 
This is to say p^ B = + p*+ B+C , p£ B = 

P& B+C +P& B+C , ^d so on. 
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FIG. 4: N = 2 . p— dependence of the normalized ratio 

p[f/p 2 = f q (p)/p 2 (top), PilVb(i-p)] = Af/W-p)} = 

[P - h(p)VW - P)\ [middle), and P <«>/(1 - P? = [1 - 2p + 
f q (p)]/(l-p) 2 (bottom). 



(17) for q = and q = 1. The AB— marginal proba- 
bilities recover Table II. For the simple particular case 
A = B = C, the Tables VII become respectively the 
Tables VIII where we have used f q (p,p) = f q (p)- 

In Fig. 6 we present Pni(p) = 2[f q (p) -p 2 } +pf q (p) for 
typical values of q and the q— dependence of p^\(2/3). 
And in Fig. 6 we exhibit the q— dependence of the cutoff 
probability p* below which the pm probability vanishes. 
Finally, in Figs. 7 we present the relevant normalized 
ratios 



TABLE VI: AB-marginal probabilities for three binary sub- 

(top) : 4 = W)- P W <( ri (26) 

The mutually independent case, specially correlated ^ ^ 

case (with p± + pf + pf > 2), and an interpolation be- , . p^i P211 P112 



(middle) 



tween them are indicated in Table VII. We verify Eq. • p 2^_ p ^ p 2 (l-p) p 2 (l-p) 
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TABLE VII: Joint probabilities for independent three bi- 
nary subsystems (top), for specially correlated three binary 
subsystems (middle), and for an interpolation between them 
(botom). 

_ 2p 2 - [M ~ Pf q (p) (07 , 
p 2 (l-p) ' [2i) 

„(«) J.i) n (i) 

(bottom) : P221 = Pl22 = ?212 

p(l — p) 2 p(l — p) 2 p(l — p) 2 

_ gjj = l-2p + f q (p) 

(1-pf (1-p) 2 • 1 ' 

III. FINAL COMMENTS AND CONCLUSIONS 

The so called "nonextensive" entropy S g (q 1) is 
in fact nonextensive (nonadditive strictly speaking) if 
we are composing subsystems A and B assumed (some- 
times explicitly, but most of the times tacitly) to be in- 
dependent. Indeed, in such case we have S q (A + B) = 
S q (A) + S q (B) + (1 - q)S q (A)S q (B), which, unless q=l, 
generically differs from S q (A) + S q (B). It is from this 
property that currently used expressions such as nonex- 
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TABLE VIII: Joint probabilities for independent three bi- 
nary subsystems (top), for specially correlated three binary 
subsystems (middle), and for an interpolation between them 
(bottom). For each cases, pi = pf — pf — p 




FIG. 5: N = 3 . Left: p— dependence of p\^ 1 for typical values 
of q. Right: q -dependence of (2/3) (p{\\(2/3) =8/27). 



tensive entropy and nonextensive statistical mechanics 
stand. However, special types of collective correlations 
do exist for which extensivity is recovered for the appro- 
priate value of q. This is to say, correlations such that 
S q (A + B) = S q (A) + S q (B). This situation has been il- 
lustrated for N = 2 and N — 3 equal binary subsystems 
in Figs. 3 and 6 respectively. 

The q = case has also been illustrated (see Q) for an 
arbitrary number N of arbitrary subsystems (with num- 
ber of states W\, W%, Wn respectively). The a priori 
total number of states is W = Y\ r=1 W r , but most of 



FIG. 6: N = 3 . Left: g-dependence of the cuttoff probabil- 
ity p* . Right: p— dependences of 3S q (l) (solid curve), and of 
S q (3) (dots) for typical values of q. The corresponding cutoff 
probabilities p* are indicated by arrows. Since numerical cal- 
culation of S q (3) becomes unstable when p is close to p* , we 
only show values of p not too close to p* . For all values of q 
such that p*{q) < 1/2, we can easily obtain the \aw-p branch 
of S q (3) just by using the permutation p with 1 — p. However, 
the situation is considerably more complex for those values of 
q such that 1/2 < p*(q) < 2/3; in this case a new branch of 
solutions is needed to cover the region 1 — p* < p < p* . We 
have not addressed this situation. 



them have zero probability. In other words, the num- 
ber of effective states (i.e., those whose probability is 

generically nonzero) is only W e ff = (J2^=i W r ) —N + l. 
Consequently, a necessary condition for this very spe- 
cial type of correlation to occur is the system to have 

W - W e ff = (]lf =1 W r ) - (E^Li W r ) + N-1 zeros in 
its table of joint probabilities. If the subsystems are all 
equal, we have W = , whereas W e ff = N(Wi -1) + 1. 
At the thermodynamic limit, it clearly is W e ff << W , 
i.e., liniN^iW^W/WiN)} = 0. If our subsystems 
were such that W r — > oo (Vr) yielding a continuum, 
they would ultimately lead to a finite Lebesgue mea- 
sure. This measure would be the hypervolume W associ- 
ated with Ndi dimensions, i.e., essentially W = JI^Li L r 
where L r is the Lebesgue measure associated with the 
d\ — dimensional r— th subsystem. In remarkable con- 
trast, W e ff would correspond to a set of zero Lebesgue 
measure, such as, for instance, a (multi)fractal whose 
Hausdorff dimension du would be smaller than Nd\ . 

Within such a scenario, it is natural to conjecture the 
following situation for q decreasing from say 1 to (see 
also When the subsystems are strictly indepen- 

dent (no correlations at all) or nearly independent (typ- 
ically short range two-body interactions within a many- 
body Hamiltonian system) , we expect an exponential N- 

dependence W ~ W e ff ~ p N (p, > 1), and the ex- 
tensive entropy to be the BG one. In contrast, if the 
subsystems have special collective correlations (typically 
long range two-body interactions within a many-body 
Hamiltonian system), we expect a power-law behavior 

W e S ~ NP > 0) << W, and the extensive entropy to 





9 = . 


pS'pV -p) = i/j>*\ 


^ q = 0.5 




q = 0.75 




q = 0,95 




/, q = 0.99 




// q ^ x - 








/ 



p ° 







/ 










9-1 \\ \ 
,/ = 0.99 — Vr 




p£/p(i -p) 1 


q - o.ys — VV — 

q - 0.75 V\ 






9 = 0-5 \ \ 






- ,=0 \ \ 






p,:;vpu-p) ! -o\ 

\, \ 





p ° 



FIG. 7: N = 3 . p— dependence of the normalized ratio 
pm/p 3 = [2(/,(p) -p 2 ) +Pf q (p)]/p 3 (top), p[^/[p 2 (l- 

p)] = pS/b 2 (i - p)] = pff 2 /b 2 (i - P)] = [2P 2 - hip) - 

P/,(P)]/[P 2 (1-P)] (middle), and pifjlpil-p) 2 ] =p { $ 2 /]p(l- 
p) 2 ]=P^ 2 /[p(l-p) 2 ]=pn^ 2 /(l-pf = (l-2p + /,(p)/(l- 
p) 2 (bottom). 



be S q with q = 1 — 1/p. Consistently, if p = 1 , then the 
extensive entropy is Sq. For all q < 1 we expect, in the 
continuum case, a zero Lebesgue measure, and a fractal 
dimension dn decreasing with decreasing q. A pictorial 
image can help understanding the conjecture. Travel- 
ling in a Brownian way — say an hypothetical "blind" 
cowboy on the back of an hypothetical "blind" horse — 
will lead to a virtually homogeneous visit of a big (rela- 
tively plane) territory, associated with a finite Lebesgue 
measure. And the result would roughly be the same in- 
dependently of the initial condition (starting point of the 
travel). This would be the typical dynamics associated 
with strong chaos (i.e., positive Lyapunov exponents), 
thus leading to the q = 1 entropy. Travelling in the way 
of a pilot of an airline company is quite different. First 
of all, he(she) will only visit the set of airports through 
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which this company operates. Given the small size of the 
airports (compared to the size of a wide territory), this 
set constitutes a set of virtually zero Lebesgue measure. 
Although statistically similar in geometrical terms, the 
result does depend on the initial conditions (the most 
important hub of the network of Japan airlines is Tokyo, 
whereas the most important of Varig is Sao Paulo, and 
the most important one of Continental airlines is Hous- 
ton) . Although no rigorous proof whatsoever is yet avail- 
able, the typical dynamics to be associated is expected to 
be that of weak chaos (i.e., basically zero Lyapunov expo- 
nents). We expect the adequate entropy to be q ^ 1. If 
we wish to recover homogeneity in the visits, we need to 
average over virtually all the possible initial conditions. 
Such an average is not needed for the q = 1 case. 

Occupancy of a phase space without strong restrictions 
makes equal probabilities of the joint system compatible 
with equal probabilities of each subsystem. This compat- 
ibility disappears if visiting of some regions of the joint 
space is strongly enhanced whereas visiting of others is 
strongly inhibited (see Figs. 4 and 7). This is clearly seen 
for the generic q — case for any value of A. It might 
be necessary to go to the asymptotic A — > oo limit in 
order to see it for say Q < q < 1. In any case, the dynam- 
ics conjectured for the q ^ 1 cases seem to be compatible 
with a recent connection Q in terms of recurrent visits in 
phase space (the q = 1 limit corresponding to a Poisson 
distribution of times between consecutive visits, during 
the time evolution of the whole system, of a given cell of 
phase space). 

Let us elaborate some more on the important connec- 
tion of q with geometry. We consider, for simplicity, 
the case of A equal subsytems Ai, A2, An, each of 
them having W\ possible microstates. The total space 
has then W = W± possible microstates that can be 
represented on a discrete A-dimensional hypercube of 
linear size W\. We shall focus however on the effec- 
tive number W e ff (W e ff < W) of microstates whose 
probability generically is not zero. A most trivial oc- 
cupation is when only the "origin" corner of the hyper- 
cube is occupied, i.e., Pxx' A i'"' An — 1- We then have 
W e ff = 1. A simple nontrivial case is when , in addition 
to that "corner" , the A "edges" of the hypercube start- 
ing from the "corner" are occupied as well (with generic 

probabilities). We then have W e ff = 1 + A(Wi - 1), 
hence p = 1 and q = 0, as already addressed. A next 
case in this series is occupancy of all the A(A — l)/2 
"faces" starting at that "corner". We then have W e ff = 
1 + N(Wi - 1) + [A(A - l)/2](Wi - l) 2 , hence p = 2 
and q = 1/2. The next case in this series is occupancy 
of all the A(A - 1)(A - 2)/6 "cubes" starting at that 

"corner" . We then have W e ff = 1 + N(W X - 1) + [A( A - 
l)/2](Wi - l) 2 + [A(A - 1)(A - 2)/6](Wi - l) 3 , hence 
p = 3 and q = 2/3. In general, if all the r-dimensional hy- 
percubes starting at that "corner" are occupied, we have 
W e ff = 1 + N(Wi - 1) + [A(A— 1)/2] (Wi - 1) 2 + [A(A - 
1)(A - 2)/6](Wi - l) 3 + ... + [A!/(A - r)!r!](Wi - l) r , 



hence p = r and q — (r — l)/r = 1 — 1/r. The last 
element of this series corresponds to fully occupy the 
unique A-dimensional hypercube. We then have W e ff = 
1 + N{Wx - 1) + [A(A - l)/2](Wi - l) 2 + ... + [A!/(A - 
r)!H](I¥i - l) r + ... + N{W X - l)^ 1 + (Wi - l) N = 
[1 + (Wi - 1)} N = hence q = 1. A different ge- 

ometry which nevertheless belongs to the q = 1/2 uni- 
versality class is the following: if, in addition to the 
"corner", all the diagonals of the A(A — l)/2 "faces" 
starting at that "corner" also are occupied, we have 
W e ff = 1 + [A(A - l)/2](Wi - 1), hence p = 2 and 
q = 1/2. Another different geometry could be the fol- 
lowing one: assume that, in addition to the "corner", 
the "edges that are (strictly or substantially) occupied 
are not all A edges, but only those following the Can- 
tor set sequence 101000101000000000101000101..., whose 
fractal dimension is d^^ = In 2/ In 3. We then have 
W H = 1 + N d H (W 1 — 1) (we are assuming that A is a 
power of 3), hence p — < 1 and q = 1— l/d^ < 0. 
A similar situation can occur for the W± states. Sup- 
pose that, in addition to the "corner", all the A "edges" 
are occupied but not fully occupied. Suppose that the 
(Wi — 1) states are fractally occupied (again the Can- 
tor, or any other sequence) with fractal dimension . 

We then have W e ff = 1 + N(Wi - l)*™, hence p = 1 
and q = 0. A quite general situation could be, in the 

(A,Wi) -» (00,00) limit, W e JJ - N d H W{ H . In 
all these illustrations, the probabilities associated with 
the W e ff occupied microstates have no particular rea- 
son for being equally probable. They could very well 
constitute a network (e.g., a scale-free network) whose 
occupancy probabilities would characterize main hubs, 
and secondary hubs, and so on (quite like the previously 
mentioned set of airports used by a given airlines com- 
pany) . In fact, it would not be really surprising if classical 
long-range-interacting many-body Hamiltonian systems 
would visit cells in phase space according to probabilities 
of this type. 

This is the typical scenario we expect for the family of 
entropies S q . It is easy to imagine that there can easily 
be even more subtle situations for which the apropriate 
(extensive) entropy would be not included in the S q fam- 
ily for any value of q. Different entropic forms would 
perhaps be then needed. But even within the S q family, 
various aspects remain to be solved, that have been only 
preliminarly addressed here. Let us mention two of them. 

First, the solutions that we have exhibited here are 
probably not unique (see the captions of Figs. 3 and 
6). Other branches of solutions could well exist. We 
have been unable, at the present stage, to find them 
all. The reader surely realizes the nontrivial mathemati- 
cal difficulty of simultaneously satisfying the impositions 
of theory of probabilities (sum of all the joint proba- 
bilities equal to unity, partial sums of the joint prob- 
abilities equal to the marginal probabilities) and those 
of extensivity of S q . The general solution seems to be 
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[l| intimately related with the recently introduced q- 
product 01 x x q y = {x x ~ q + y x ~ q — l) 1 /^ 1 -?^ which has, 
among others, the following properties: (i) x x i y = xy, 
(ii) hi q (x x q y) = ln 9 x + ln q y (whereas ln q (xy) — 
\n q x + ln q y + (1 - q) In, x In, y) ; (iii) l/(ic x q y) = 
{I fx) x 2 - q (1/y); (iv) x x q (y x q z) = (x x q y) X q z = 
xx q yx q z = {x l - q +y 1 - q +z 1 -i-2) 1 /( 1 -^- ) (v) xx q l = x. 
This interesting structure probably is one of the ingredi- 
ents, but there are surely others to be considered con- 
comitantly, specifically those related to the constraints 
imposed by theory of probabilities. 

Second, to illustrate an important point let us rewrite 
Eq. (22) as follows: 

2p K + 2(1 -pT - (f K ) K - 2(p - /„)" - (1 - 2p + f.T = 1 • 

(29) 

This relation means that it has been possible to find a 
function f K (p) which satisfies the impositions of theory of 
probabilities. Now, if we wish this solution to correspond 
to the extensivity of S q , we just impose k = q. If we 
wish instead to impose the extensivity of S2- q we identify 
k = 2 — q. In this case, we have solutions corresponding 
to 1 < q < 2. We can even impose, if we wish, the exten- 
sivity of S K (g\ , where n(q) is virtually any (increasing or 
decreasing) monotonic function of q satisfying k(1) = 1. 
This freedom might play a relevant role when a thermody- 
namical (or thermodynamical-like) connection is seeked. 
Indeed, most of the systems which seem to obey nonex- 
tensive statistical mechanics exhibit a (quasi) stationary 
state whose entropic index is q > 1. This point needs 
further analysis in order to unambiguosly establish the 
identification between k and q which is thermodynami- 
cally adequate. It is interesting at this stage to recall a 
recent discussion by Robledo [6j on a nonthermodynam- 
ical system, which has nevertheless some analogy with 
the present situation. Basically, the so called Mori's q- 
transitions for say the usual logistic map occur at both 
.2445... and 2 - 0.2445... 0. 

A few words on terminology to conclude. We have 
seen that (under specially correlated composition of sub- 
systems) S q can be strictly additive (i.e., S q (A + B) = 
S q (A) + S q (B)) for a variety of values of the entropic in- 
dex q. It has nevertheless become current denomination 
to refer to the q = 1 universality class as extensive or 
normal systems, and to the g ^ 1 universality classes as 
nonextensive or anomalous systems. This use has origi- 
nated from the fact that, in the early times of the the- 
ory, the focus was explicitly or tacitly put onto indepen- 
dent subsystems. For this simple composition law, and 
only then, Sbg * s strictly additive (i.e., Sbg(A + B) = 
Sbg(A) + Sbg(B)), whereas S q (q ^ 1) is not (i.e., 
S q (A + B) = S q (A) + S q {B) + (1 - q)S q {A)S q (B)). Al- 
though slightly misleading from the entropy standpoint, 
the current notational distinction extensive versus nonex- 



tensive is instead perfectly natural from the energy stand- 
point of Hamiltonian systems. Indeed, suppose we have a 
d-dimensional classical system with attractive two-body 
interactions whose potential energy decays with (dimen- 
sionless) distance r as —l/r a (a > 0). Let us further 
assume for simplicity that a strong repulsion exist at 
the r — > limit (therefore no nonintegrable singulari- 
ties exist at short distances). The Lennard- Jones gas 
would be (a, d) = (6,3); Newtonian gravitation would 
be (a, d) = (1,3) (if we take into account the fact that 
at very short distances, important repulsive quantum ef- 
fects are expected which would avoid the mathematical 
problems tied to the — 1/r singularity). The total poten- 
tial energy at the ground state is expected to be U(N) oc 

-N if 1 '" drr d - x r- a oc N[N 1 ~ a / d -l]/(l-a/d), where 
we have assumed for simplicity that the N elements of the 
system are roughly homogeneously distributed in space, 
and where the dimensionless distance r has been chosen 
to be unity at the short distance effective cutoff. We im- 
mediately see that the energy is extensive if a/d > 1, 
whereas it is nonextensive if < a/d < 1. It is long 
known (see, for instance, @) that, for the a/d > 1 
systems, BG statistical mechanics is perfectly adequate. 
More over, for them the t — > oo and the N — > oo limits 
commute, thus always leading to thermal equilibrium. 
On the other hand, plethoric evidence now exists that, 
in remarkable variance, for the < a/d < 1 systems, 
the J -> oo and the N — > oo limits do not commute, 
the physically interesting states for large systems being 
the (quasi)stationary or metastable ones corresponding 
to taking first the N — > oo limit and only afterwards the 
t — > oo limit. For such anomalous states, the inadequacy 
of BG statistical mechanics is notorious when we use no 
other dynamics than the natural one (Newton's law if the 
system is classical) 9]. For instance, the distribution of 
velocities is seemingly not Maxwellian, and there is aging, 
a phenomenon absolutely incompatible with the transla- 
tional invariance expected for thermal equilibrium. A 
transparent proof that nonextensive statistical mechan- 
ics (with a 9 / 1 entropy S q ) is in place is still lack- 
ing, but this could be the case. Indeed, vanishing Lya- 
punov exponents have been exhibited, as well as the spe- 
cific anomalous diffusion associated with the nonlinear 
Fokker-Planck equation, and a variety of q— exponential 
behaviors Q. Work is in progress and further contribu- 
tions are welcome. 
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